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Abstract 

We investigate nuclear spin-lattice relaxation data in the normal state of optimally doped 
YBa2Cri307 by analyzing the contributions to the relaxation rate of the copper, planar oxygen 
and yttrium along the directions perpendicular to the applied field. In this new picture there is 
no contrasting temperature dependence of the copper and oxygen relaxation. We use the model 
of fluctuating fields to express the rates in terms of hyperfine interaction energies and an effec- 
tive correlation time T e g characterizing the dynamics of the spin fluid. The former contain the 
effects of the antiferromagnetic static spin correlations, which cause the hyperfine field constants 
to be added coherently at low temperature and incoherently at high temperature. The degree of 
coherency is therefore controlled by the spin-spin correlations. The model is used to determine 
the temperature-dependent correlation lengths. The temperature-dependent effective correlation 
time is found to be made up of a linear and a constant contribution that can be related to scatter- 
ing and spin fluctuations of localized moments respectively. The extrapolation of our calculation 
at higher temperature fits the data also very well at those temperatures. The underdoped com- 
pounds YBa2Cu306.63 and YBa2Cu40s are studied in the limit of the data available with some 
success by modifying the effective correlation time with a gap parameter. The copper data of the 
La2_ x Sr a; Cu04 series are then discussed in terms of the interplay between the two contributions 

to T eff . 

PACS numbers: 74.25.-q,74.25.Ha,76.60.-k 
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I. INTRODUCTION 



The nuclear resonance techniques, nuclear magnetic resonance (NMR) and nuclear 
quadrupole resonance (NQR), are powerful probes for investigating the microscopic mag- 
netic properties of cuprates that exhibit high-temperature superconductivity. An eminent 
advantage of these methods is based on their highly local nature which allows one to get 
information about the distinct chemical species in the materials and make selective measure- 
ments of different crystallographic sites (for reviews, see Refs. Qj^j^P). Magnetic hyperfme 
interactions couple the nuclear spins to the electron system and it is essential that they are 
known as accurately as possible in order to allow a correct interpretation of the properties 
of the electron liquid in terms of measured NMR or NQR data. 

One particular aspect that absorbed the attention of the NMR community is that the exper- 
iments in YBa 2 Cu 3 7 , YBa 2 Cu30 6 6 3, and in Lai. 85 Sro.i5Cu0 4 (Refs. 

seem to show a dramatic contrast in nuclear-spin lattice relaxation rate behavior between 
copper and oxygen sites in the CuCVplane, although these sites lie less than 2 A apart. 
It was concluded that the relaxation at the copper site exhibits strong antiferromagnetic 
(AFM) enhancement effects, whereas that at the planar oxygen site is weakly enhanced 
with strikingly different temperature dependence. 

The nuclear spin-lattice relaxation rate ^Tj^, 1 for a nuclear species k is the rate at which 
the magnetization relaxes to its equilibrium value in the external magnetic field applied in 
direction a. The relaxation of the nuclei under consideration in the cuprates is caused by 
two or more fluctuating hyperfme fields that originate from magnetic moments localized 
near the coppers. Since the squares of these fields come into play, one of the first tasks 
when interpreting spin-lattice relaxation data is to determine whether the hyperfme fields 
should be added coherently or incoherently at the nucleus. While Mila and Rice^ added 
them incoherently, Monien, Pines and Slichter— considered both extreme cases, and from 
the analysis of the copper data in YBa2Cus07 concluded that within a one-component model 
the fields should be added coherently. The question of coherency was put aside when Millis, 
Monien and Pines (MMP) gave a quantitative and complete phenomenological description 
of the relevant measurements by putting forward a models where the nuclear spin-lattice 
relaxation rate, via the fluctuation-dissipation theorem, is expressed in terms of the low- 
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frequency limit of the imaginary part of the spin susceptibility x 

<? 

The form factors k F a (q) depend on the geometrical arrangements of the nuclei and the 
localized electronic spins. Under the form Eq. the question of the degree of coherency 
is delegated to the choice of a form for the susceptibility The "basic" idea behind the 
MMP model was to account for strong AFM correlations which exist in the cuprates even 
in the overdoped regime. The MMP model therefore postulated a spin susceptibility which 
is strongly peaked at the AFM wave- vector Q = (jr, tt). In this way, the seemingly different 
relaxation behavior of copper and oxygen could be understood and almost all NMR and 
NQR relaxation measurements in cuprates have been analyzed using the MMP approach. 
In a later development of the modeU^, the susceptibility u>) was split into two parts, 
X — Xaf + Xfl- The first term, xaf represents the anomalous contribution to the spin 
system and is peaked at or near Q. The second term, xfl, is a parameterized form of 
the normal Fermi liquid contribution. For copper, the contributions from the Fermi liquid 
part are much smaller than those from \af but they dominate the relaxation of oxygen 
and yttrium nuclei. In the parameterization for \ introduced by MMP, xaf is strongly 
peaked at Q and the hyperfine fields at the copper are added essentially coherently. On 
the other hand, for any ^-independent xfl, the contributions of these fields are strictly 
incoherent, and the fields at the oxygen are therefore added incoherently in the MMP theory. 

The goal of our work is to get the information on the degree of coherency directly 
from the data. We therefore intentionally avoid the use of Eq. since the degree of 

coherency of the hyperfine fields at a nucleus depends on details of the susceptibility which 
are difficult to model. Moreover, the possibility that two or more processes contribute to 
the spin relaxation might make it difficult to disentangle these different contributions. We 
therefore stay in the direct space since there are only a few points in the lattice which are 
relevant for NMR. In this case the coherency is related to the spin-spin correlations, which 
we take as a parameter that we deduce from the experiments. The spin-spin correlations 
are found to be temperature-dependent, and therefore the degree of coherency varies with 
the temperature. As expected, the coherency is reduced when the temperature is increased. 
We note that this is also what MMP finds in the case of the copper nuclei: the peak in the 
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susceptibility (as well as the sometimes forgotten cut-off^) gets broader as the correlations 
decrease. However, the approach taken here reveals other advantages: in contrast to MMP, 
the relaxation of both the copper and the oxygen (as well as the yttrium) is caused by 
the same temperature-dependent relaxation mechanism of the spin liquid, which we will 
characterize by an effective correlation time. Also, the different temperature behavior of 
the oxygen and copper nuclear relaxation rates can then be explained naturally from the 
particular values of the hyperfine field constants. 

The premises of our different approach to the analysis of spin-lattice relaxation data 
in the normal state of cuprates are the following. We go back to the simple model of 
fluctuating fields 18 that has been applied by Pennington et al. 19 to analyze their spin-lattice 
relaxation rate data for copper in YP^CuaCV This model is particularly appropriate for 
anisotropic substances, as it is the case of the cuprates. We retain, however, the AFM 
correlations which are an essential feature of the concept of the nearly AFM Fermi liquid. 
The normalized AFM spin-spin correlations between adjacent coppers are a key quantity, 
since they determine to what extent the hyperfine fields are added coherently. 
We adopt the usual form for the spin Hamiltonian for copper as proposed by Mila and 
Ricei^, whereas that for the oxygen is determined by transferred contributions from the two 
nearest neighbor copper ions (Shastry^). We note that quantum-chemical calculations^^ 
have shown that contributions to the oxygen hyperfine interaction arising from further 
distant copper ions are marginal and that an introduction of a substantial transferred field 
from next nearest neighbor Cu is not justified. 

We will assume that all data obtained in the normal state of the cuprates can be attributed 
to purely magnetic relaxation although there are indications^ that some observed phenom- 
ena hint for an additional influence of charge fluctuations. 

We confine our analysis to measurements of the planar copper and oxygen and the yttrium 
nuclei in the normal state. We neglect orthorhombicity and assume a quadratic CuCVplane 
with a lattice constant of unity. For the planar oxygen we distinguish between the direction 
a parallel to the Cu-O-Cu bond, and b, perpendicular to the bond. 

Our work is structured as follows. We introduce in Sec. |H] a new representation of the 
relaxation data which is appropriate to their analysis in anisotropic materials. Sec. 11111 
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presents the model and assumptions made. We then apply our treatment to YBa2Cus07 
in Sec. IIVI and confront our results with a range of experiments, in particular those made 
at high temperature. The significance of the model parameters at low temperature are 
discussed. Due to lack of data, the underdoped materials YBa2Cu 3 6 . 6 3 and YBa 2 Cu 4 08 
are less extensively treated in Sec. |V| We turn to the La 2 _ x Sr x Cu04 series in Sec. I VII and 
discuss the temperature and doping dependence in view of our model, in particular in the 
high temperature limit. We present a summary and some conclusions in Sec. IVHI Some 
special considerations are collected in the Appendices. 



II. REPRESENTATION OF DATA FOR ANISOTROPIC MATERIALS 

In this section we plot the spin-lattice relaxation data in a way which suggests that the 
same relaxation's mechanism is at work for all the nuclei under consideration. 
The magnetic relaxation process into equilibrium is caused by fluctuating effective magnetic 
fields along the two orthogonal axes (3 and 7 perpendicular to the direction a of the applied 
field, which we write as 

k T£ = k U p + % (2) 

where (a, (3, 7) = (x,y,z). The quantities k U a describe then the contribution to fc T 1 ^ 1 and 
k T{ J 1 caused by fluctuating fields in the crystallographic direction a. From (J2J) the k U a are 
therefore given by 

k U a = \ P'TfJ + %J + (3) 

These transformed rates k U a , which in the following will just be called rates, are not directly 
accessible by experiment except for particular symmetries (e.g., 63 U a b = 63 T 1 ~ 1 /2), but in 
general can be obtained if a complete set of data measured with the applied field along all 
three crystallographic axes is available. At variance with the common use to represent NMR 
spin-lattice relaxation rates in (TiT) -1 versus T plots, we prefer here to study the rates U a . 
This change in representation is trivial. It was however pointed out previously by one of 
us^ that it is much more instructive to investigate ratios between k U a for a in different 
directions than ratios of the corresponding rates k T^ a . 

The practice to represent and analyze (TiT) -1 data grew out from applications in liquids 
where the rates are isotropic and (TiT) -1 is temperature independent in simple metals. For 
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layered cuprates, however, which are anisotropic materials, the analysis of NMR data in 
terms of k U a has distinct advantages over that in terms of k T± a . 

For cuprates, the most complete set of NMR and NQR data with respect to different 
nuclei and directions of external field relative to the crystallographic axes is available for 
optimally doped YBa 2 Cu307 in the normal state between 100 K and room temperature. We 
start, therefore, our analysis in this temperature range and use the copper data for 63 T 1 ~ 1 
from Hammel et air and 63 7\~^ from Walstedt et a/.—. The oxygen data 17 T 1 ~ 1 were taken 
from Martindale et alM and the yttrium data 89 from Takigawa et alM-. 

These sets of data for the rates k T{^, interpolated at the same temperature points, allow 
the transformation into our rates k U a according to 0. Neglecting any error analysis, the 
results for the three nuclei are shown in Fig. ^ 

In these representations of the data there is nothing to see of a drastic contrast in the 
temperature dependence of the relaxation rates of copper and oxygen nuclei in the same 
CuC>2 plane which has, as mentioned in the Introduction, intrigued the NMR community. 
All relaxation rates grow with increasing temperature as is expected for fluctuations. The 
differences in the magnitude for Cu, O, and Y are due to different strengths of the hyperfine 
interaction energies which allow a scaling of all data as is shown in Appendix [XJ We stress 
that so far no model is used, and nothing else was done other than adding and subtracting 
the data. 

The relaxation rate data, transformed now into a representation appropriate for layered 
structures, suggest therefore that all nuclei under consideration relax in a similar fashion 
by the same mechanism of the spin liquid and that the strength of this relaxation initially 
increases linearly with temperature (note that the additional line drawn in Fig. |1 (a)] depicts 
the relaxation measured in copper metal). The question of why the relaxation rates between 
the copper and the oxygen are nonetheless different will be investigated in IIV Fl In the 
following section we outline a robust and unbiased model that will allow us to deduce more 
details about the origin and the source of the temperature dependence of the spin-lattice 
relaxation. 
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FIG. 1: Relaxation rates k U a from measurements in YBa2Cu30"7. The symbols denote interpola- 
tions of T^" 1 dat a 8 i 25 i 26 i 27 transformed according to (JHJ), in the direction a (diamonds), b (squares), 
ab (circles), c (triangles). Dashed lines are guides for the eyes. The full straight line in (a) denotes 
the values for metallic copper. 



III. MODEL AND ASSUMPTIONS 



The quantities U a have been introduced as the contributions from fluctuating local 
effective magnetic fields H a along the direction a. In this section we adopt the calcula- 
tion of relaxation rates in terms of fluctuating fields described in Ref. [l8| to determine 



an expression for k U a in the present case. In particular in the cuprates the hyperfine 
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FIG. 2: Schematic picture of the AFM spin correlations between electron moments at the 
copper sites 0, 1,.. 4. 



interaction energies depend on the static AFM spin correlations. We find that k U a , 
for any nuclei k under consideration, can be written as a static term containing the spin 
correlations and hyperfine field constants, and another which is the effective correlation time. 

Let us consider first an oxygen nucleus with spin 11 1 and gyromagnetic ratio 17 7, for 
which the hyperfine interaction Hamiltonian is determined by 

17 H(t) = -h 17 7 17 H(t) • 17 1 . (4) 

The field operator 17 H(t) for an oxygen situated between site and 1 (Fig. |2J is assumed to 
originate from magnetic moments with spin S = 1/2 localized on those two adjacent nearest 
neighbor (NN) copper ions with spin components S® and respectively (Fig. |2J). The field 
operator components are given by 

17 H a (t) = -h le c a [S° a (t) + S 1 a (t)} (5) 

where c Q is the diagonal element of the hyperfine tensor in direction a given in units of spin 
densities (a^ 3 ). 

Provided that the fluctuations in different field directions are independent, the components 
of the auto-correlation function of 11 H(t) are 

( 17 H a (t) 17 tf Q (0)> = ff^cl 2 {S%(t)l%(P) + SS(t)S?(0))- (6) 

Assuming an exponential decay of correlations in time, the expression in brackets can be 
written as 

(SZ(t)S%(0) + S%(t)SnO)) = \{l + e-NA*. (7) 
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r e ff is an effective correlation time that acts as a time-scale for the fluctuations^. Its 
temperature dependence will give us an indication as to what type of processes come into 
play in the relaxation. is the normalized "static" NN spin correlation 

7^ = 4(^(0)^(0)). (8) 

The correlations of two NN moments is thus separated into a spatial AFM correlation 
that will determine the degree of coherency and temporal correlation that characterizes the 
dynamics of the electronic spin fluid system which exchanges energy with the nuclei. 
We make the assumption that the correlation time r ej gr is isotropic and in the following that 
it is also independent of the spatial separation of the correlated spins. However may be 
different for in-plane and out-of-plane components. Following Slichter 18 , we find that k U a is 
obtained as 

1 f°° 

kUa= 2Vj ( kH ^) k H a (0))e luJkt dt (9) 

where Uk is the Larmor frequency. The expression for the contribution 17 U a to the relaxation 
rate then is given by 

U Ua = J^(l + Kl)Teff (10) 

when \ukT e ff\ «C 1, with C a = h 17 ~f h^ e c a . 

The hyperfine field at the planar copper nucleus is determined by an on-site contribution 
A a from the copper ion with spin component Sq and transferred contributions B a from the 
four NN ions with spin components S", where j = 1, 2, 3, 4 (Fig. |2J). We note that ab-initio 
calculations of the hyperfine interactions^ 1 ^ yield besides the isotropic transferred field B 
also a transferred dipolar field Bdi P , a which, for simplicity, will be ignored in the following. 
The corresponding equation for 63 U a then reads 

63 U a = ^ (A 2 a + AB 2 + 8A a BK^ + 8B 2 K? 2 + AB 2 K^) r eff (11) 

where K" 2 and K® 3 are the normalized spin correlations between two copper ions which are 
y/2 and 2 lattice units apart respectively (Fig. |2J). 

In the YBaCuO compounds, the yttrium is located between two adjacent CuC>2 planes 
and has four copper NN in each. Here, the situation is not as clear as for the Cu and 
O hyperfine interactions since there may be interplane spin correlations between copper 
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moments and a direct dipolar coupling of the same order magnitude as the transferred 
fields. We ignore these complications and use the simplest form which leads to 

89Ua = \3 (1 + 2Kl + K?2) Teff - (12) 

It is convenient to express the relaxation rates as 

k U a (T)= k V a (T)r eff (T) (13) 

with 



17 V a (T) = i ^2C2[l + AJ 1 (T)] (14a) 
63 V a (T) = ± [A 2 a + AB 2 + 8A a BK^(T) + 8B 2 K? 2 (T) + 4B 2 K? 3 (T)} (14b) 
89 V a (T) = ±- 2 8D 2 a [l + 2Kg 1 (T) + K? 2 (T)}. (14c) 

This factorization emphasizes the different temperature dependencies that determine 
k U a (T). The k V a (T), which apart from the factor h 2 are the static hyperfine energies squared, 
vary with temperature due to changes of the static spin correlations K°-(T), whereas the 
effective correlation time r e ff(T) reflects the changes in the dynamics. 
We denote the limiting values when all correlations are zero by k V® 

^ = 4^K + 4 ^ 2 ) (!5a) 

17 V;° = ^C 2 a (15b) 

8 V a ° = ^8D 2 a (15c) 
and by k V[ c for full AFM correlations: 



a 

17 



6 M C = T^(^" 45 ) 2 ( 16a ) 

V£ = (16b) 

m V[ c = 0. (16c) 

k V a (T) accounts for the transition from the fully correlated situation, where the hyperfine 
fields are added coherently, towards the completely uncorrelated regime given by fc V Q °, 
where the fields are added incoherently. 
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We would like to point out that the model we use is an extension of a well established 
approach. It has been applied, e.g. by Monien, Pines and Slichter— to analyze 63 Xf 1 data 
in the limiting cases 63 V and 63 V* c . The new feature is to adopt the existence of AFM 
correlations in the cuprates that are static with respect to typical NMR times and vary 
with the temperature, and this also has been done before, as it is detailed in section II V Gl 

In order to reduce the number of parameters we assume that the static spin correlations 
are antiferromagnetic and that K° 2 and K® 3 depend directly on the value of Kq 1} according 
to 

Kt 2 =\K^ and =\ | 2 . (17) 

This particular choice of an exponential decay with length X a , defined so that 

i^ 1 = -exp[-l/A a ], (18) 

is guided by results of solutions of the planar anisotropic Heisenberg model which were ob- 
tained by direct diagonalization of the Hamiltonian for small systems^ and which suggested 
Kf 2 ~| | 1-5 and anisotropic correlations. 

These assumptions certainly influence the actual values for A Q which will result from the 
analysis of the data, but they provide a reasonable starting point for the analysis and may 
be easily modified. 

IV. RESULTS FOR YBA 2 CU 3 7 

We present in this section the analysis of the optimally doped YBa2Cus07 data in the 
framework of our model. We first outline the fitting procedure, which allows us to determine 
the correlation lengths A a b and A c , and the hyperfine field constants. We can then determine 
the effective correlation time T e g and parameterize it in order to identify the underlying 
mechanisms of the spin fluid. We discuss in particular the low temperature limit, which 
exhibits a Fermi-liquid character for all three nuclei (Cu, O and Y). The model predictions 
at high temperature are then compared with experiments. Then, by looking at the extremal 
values of k V a , we will discuss why the measured relaxation rates of the copper and oxygen 
have different temperature dependence. Finally, we apply our model to spin-spin relaxation 
measurements. 
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A. Fitting procedures 



We define the following independent ratios rj 



63 u c _™U ab _™U ab ^U b 89 f/ c 17 u c 

ri ~^U c T2 -^lT a r3 ~^Lh U ~W e r5 "^ r6 -^U c (19) 
and denote with r exv the ratios obtained from the experimental values for k U a as calculated 

from (J3J) and plotted in Fig. ^ We also form the ratios r™ ,od according to the model (Eqs. 
110111112]) . for which the correlation times r e g cancel (r™ od = 63 V c / l7 V c etc). 
We define the following function to minimize: 



x 2 



where n r (6 in the case of YBa2Cus07) is the number of ratios available and j runs through 
the temperature points. The function measures the normalized difference between the cal- 
culated ratios rf- od and the experimental points r^ xp at each temperature. We then minimize 
(120)1 in order to obtain the best local solution for the whole set of parameters (hyperfine 
interaction energies and values for the correlation lengths). As input parameters for the 
hyperfine interaction energies we used the values determined by Hochnep^. The X ab and 
A c obtained for each temperature point are illustrated in Fig. El and the resulting values 
for the hyperfine interaction energies^ 8 - are given in Table |U It is seen in Fig. E] that A c is 
consistently somewhat larger than \ ab (A c ~ l-2A a b) which could be related to our neglect of 
anisotropies in the g-factors. From the logarithmic representation in the inset of Fig. El we 
see that the temperature dependence obtained for the X a can be parameterized by a sum 
of two exponential functions. This parameterization has been only used to extrapolate at 
higher and lower temperatures (dotted lines Fig. |4(a)| ). We would like to emphasize that 
we are not at this stage looking into interpreting the temperature dependence of A a . We 
merely wish the A Q to be optimally fitted (regardless of the significance of the parameters), 
and this particularly at high temperature. 

The NN spin correlation functions are plotted in Fig. |4(b)| The correlation Kq\ is 
about —0.4 at T= 100 K and drops to —0.2 at room temperature. The extrapolation to 
higher temperatures shows that even at T = 600 K small AFM correlations (~ —0.04) exist. 
It is clear that these specific values depend on the particular choice (Eq. fTT|) of exponential 
dependence of the correlations with the distance. 
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FIG. 3: Correlation lengths A a obtained from the fit in YE^CusOy, in the direction ab (circles) 
and c (triangles). Inset: the same on a logarithmic scale, fitted with exponential functions. 
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FIG. 4: Correlations and correlation lengths from the fit (lines) in the direction ab (circles) and 
c(triangles). The dotted lines are extrapolations with exponential functions. 



We note that the hyperfine interaction energies in Table H] have been only determined 
from relaxation data without recourse to static measurements except that some reasonable 
initial guess had to be assumed. It is astonishing that the resulting values are very close to 
those which have been compiled by Nandor et alM and which are included in Table 

For further reference we collect in Table HT1 the values obtained for k V® calculated from 
(|15|) using the hyperfine interaction energies given in Table |l] 
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63 A c 


63 A ab 




17 c a 


17 c b 


l7 c c 


m D ab 


89 D C 


Fit 


-1.68 


0.168 


0.438 


0.259 


0.173 


0.196 


-0.00280 


-0.00349 


Ref^ 


-1.6 


0.29 


0.4 


0.25 


0.13 


0.156 


-0.0048 


-0.0048 


"Rcf. 22 


-1.72 ±0.42 


0.21 ±0.11 


0.29 0.05 


0.289 


0.140 


0.151 







"Cluster calculations. The error intervals are due to uncertainties in the spin-orbit interaction. 



TABLE I: Magnetic hyperfine interaction energies in units of 10~ 6 eV 



17y0 
v a 


17 v b ° 


17y0 
' c 


63t/0 
V ab 


63 yO 


89t/0 
v ab 


89 yO 


77.7 


34.6 


44.2 


459 


2080 


0.0361 


0.0562 



TABLE II: k V° in units of 10 15 s" 2 



B. Consistency check 

In the previous subsection we have formed the ratios rf ,od in order to get rid of the effective 
correlation time T e jj in the model. Having now the fitted values for the hyperfine interaction 
parameters and for Kq X (T) and thus calculated the values k V a (Eqs. flljl . we get back to the 
experimental spin-lattice relaxation rates h U a . From (|13p we can obtain the resulting seven 
experimental values k T e jj ta , which, according to the model assumptions, should all be equal. 
These seven values for T e g are plotted in Fig. EJin function of the temperature. 

The consistency is surprisingly good. We therefore calculated the mean values (circles 
in Fig. neglecting any weighting according to the estimated precision of the data for the 
various nuclei). 

The temperature dependence of the average r e g is well fit by a function of the form 

t^=t^ 1 + t 2 1 (21a) 
with Ti = aT and r 2 = constant. (21b) 

We obtain a = 7.0 x 1CT 18 s/K and the temperature independent r 2 = 3.0 x 10~ 15 s. 
The result of the fit is shown in Fig. |H1 together with extrapolations to lower and higher 
temperatures. 
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(thin dashed), 63 r a f, (thick dash-dotted), 63 r c (thick dashed), 89 r a f, (thick dash-dot-dot-dashed), 
89 r c (thick dot-dash-dash-dotted). 
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FIG. 6: Temperature dependence of the averaged relaxation time T e g. The circles are the calcu- 
lated average, and the dashed lines represent the upper and lower bound. The dotted line is the 
fit as explained in the text. 

C. The "basic" relaxation mechanism 

We would like to point out that our analysis reveals something quite different from any 
previous analysis. The fact that all the r a effectively fall onto the same line demonstrates 
that the relaxation of all the three nuclei under consideration is governed by the same 
mechanism. 
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In our picture, there exists what we will call a "basic" electronic relaxation mechanism, 
caused by the spin fluid, that affects the localized moments and exchanges energy with 
the nuclei. This mechanism is characterized by the short effective correlation time r e g. In 
addition, on a long time scale, the spin correlations between these moments vary with the 
temperature. The observed nuclear spin-lattice relaxation rate k U a reflects the temperature 
dependence of both the "basic" electronic relaxation and the change in spin correlations. 
In our model, we can disentangle these two contributions by introducing a "basic" nuclear 
relaxation rate k U a (T), defined by 

k U a (T) = k U a (T) = k V° r eff (T). (22) 

The temperature dependence of the "basic" nuclear rate k U a (T) is therefore that of the "ba- 
sic" electronic relaxation mechanism by definition. Important information on the electronic 
system is of course also drawn from the temperature dependence of the spin correlation 
K^yiT) (and hence of k V a (T)) and will be discussed in IIV Fl 

From the consistency check (|IV Bj) . we therefore conclude that the temperature dependence 
of the "basic" relaxation rate 

(i) is the same for Cu, O and Y 

(ii) is linear in T at low temperature. 

D. Two "basic" electronic relaxation mechanisms 

There are four points which would now require a discussion. The correlation time T\ 
that dominates the rates k U a at low temperatures, the correlation time Ti that dominates 
at high temperature, the crossover between the two regimes, and the particular form of T e g 
(Eq. I21a|) . where the rates are combined. There are various possible explanations of the 
crossover from the initially linear behavior to a saturation at high temperature. We confine 
ourselves to some "basic" and simple arguments in order to avoid an over-interpretation at 
this stage. 

We address first the low temperature regime and consider only the contribution T\ = aT. 
This corresponds to the behavior of T-f 1 in a metallic system. Pines and Slichter— have 
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discussed magnetic relaxation by fluctuating fields in terms of a random walk approach 
and considered also the nuclear relaxation by conduction electrons in metal. Adopting their 
argumentation to the present case, the time T\ is interpreted as T\ = p(T)rd we ii- This dwelling 
time r d W eii = fr/Ep is roughly the time a conduction electron spends on a given atom, and 
p(T) denotes the probability that during T dwe ii a nuclear spin flip occurs. For a degenerate 
Fermi gas, this probability is p(T) = T /Tp where Tp is the Fermi temperature. With these 
arguments we can express our parameter a as 



k B T F 2 v ' 

which yields for the Fermi temperature Tp = 1050 K. Moreover, we find then that the 
residence time on an atom is r dwell — 7.3 x 10~ 15 s. The rather low valued of Tp indicates 
that the degeneracy is lifted. Instead of a well defined Fermi edge, the distribution is 
smeared out and the temperature dependence of the chemical potential becomes important, 
preventing a quantitative analysis without further information. 

Adapting these arguments to the present case, it means that at low temperature all 
the spin-lattice relaxation rates (and that also for Cu) in YBa2Cu307 can be explained by 
scattering from quasiparticles within a kind of Fermi liquid model. The influence of the 
AFM spin correlations is however manifest in the temperature-dependent interaction k V a 
whose properties will be investigated in subsection IIVF1 



At high temperatures the data can no longer be explained by a Fermi-liquid behavior 
since the influence of the contribution r 2 grows. According to the fit values, we have T\ = r 2 
at 420 K. We defer a discussion of r 2 and the crossover T\ to r 2 to later sections, until we 
have looked at more NMR/NQR data in other cuprates. 

As concerns the ansatz ()21aj) for the effective correlation time we note that the interactions 
between the nuclei and the electronic system are the same whatever the "basic" relaxation 
mechanisms are. They are determined by the hyperfine energies k V a (T) which, however, 
change from strong AFM correlations at 100 K to weak AFM correlations at around 400 
K. Whether we really have two different "basic" relaxation mechanisms at work or whether 
they are two different manifestations of the same mechanism is an open question. 
A possible explanation for the special form of ()21a|) is as follows. Adding two rates T~ hort and 
T kmg means that the phase correlations between the processes associated with the longer time 
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are destroyed by the impacts of those associated with the shorter time, since immediately 
after an impact of a process with time T\ ong another impact of a process with the short time 
7 'short occurs. In our case, at low temperature r short = T\ and T[ ong = r 2 , and vice- versa at 
high temperature. In addition, we also note that r e g was originally introduced in Eq. (J7J) in 
order to express the time dependence of {Sq (fyS^ (0)). If the time evolution is now governed 
by two processes so that the Hamiltonian is = S)i + Sj 2 , we get instead of Eq. (J7J) that 
(Sq (t)S'o (0)) oc exp (— \t\/rx — \t\/T 2 ), provided Sji and S) 2 commute. We see that r e g would 
thus be given by Eq. (j21a|) . 

E. Comparison of model predictions and measurements at higher temperatures 

In YBa2Cu3C>7 several spin-lattice relaxation rate measurements have been performed 
at temperatures well above the temperature interval which we used for fitting the model 
parameters. So far we have not made use of these data since they do not allow us to 
extract the individual contributions U a . Now, however, it is of interest to compare the 
extrapolated theoretical predictions with these data. We return to the customary ( fc T lQ T) 1 
representation and show in Fig. |7(a)| a plot of ( 63 Ti a T) 1 versus T. The high-temperature 
data of Barrett et alM- are denoted by triangles down, the original data^^ used for the fit 
below room temperatures by full circles and full triangles and the model fit by the line. In 
Fig. |7(b)| we depict the temperature dependence of ( 17 Ti Q T) 1 and include data points from 
Nandor et al. 29 who also reported values for ( 89 T lc T) 1 which are shown in Fig. |7(c)[ 

The predictions of the model obtained by extrapolation to higher temperatures are in 
excellent agreement with the measured relaxation rates of the copper. The downward trend 
of the relaxation rate of the oxygen is also well reproduced. Less good agreement is achieved 
in the case of the yttrium. The experimental rate seems to be more or less constant from 200 
K onwards, whereas we predict a slowly decreasing rate. This might be due to the simplistic 
model (Eq. IT2*j) we use, neglecting interplane spin correlations and dipolar couplings. We 
also note that there is some disagreement between the data of Nandor et alM and those of 
Takigawa et ali^ already in the temperature range below 300 K. 

On the whole however the model seems to explain the general trends of the data well. A 
close inspection of the oxygen data, however, shows that the observed increase in the ratio 
of 17 T la 1 / 17 T 1 ^ 1 with decreasing temperature reported by Martindale 26 is not reproduced. 
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FIG. 7: Full symbols: original experimental data. Lines: fit as explained in the text. Hashed 
symbols: other data (see text). Direction a (diamonds), b (squares), ab (circles) and c (triangles). 

The increase of the ratio 17 T 1 ~ 1 / 89 T 1 ~ 1 for decreasing temperature observed by Takigawa et 
a/.— is qualitatively reproduced but lacks quantitative agreement. The ratio 63 T 1 ^, 1 / 17 T 1 ~ 1 is 
however in very good agreement with our predictions and will be discussed in Section IV Al 

F. Temperature dependence of interaction energies and Korringa relation 

It is instructive to discuss the calculated values k V a (T) shown in Fig. |HJ The bars on 
the left in each figures indicate the extremal values k V® (Eqs. ITrj|) if the system was fully 
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FIG. 8: Temperature dependence of k V a in the direction a (diamonds), b (squares), ab (circles) and 
c (triangles), which reflects the change in the degree of coherency. The bars at high temperatures 
indicate the values k V® with vanishing correlations, those at T = the fully correlated values k Vi 
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antiferromagnetic. The bars on the right show the other extreme k V® (Eqs. IT5j) . for a system 
with no correlations left. In the case of the oxygen (Fig. |8(b)D , 17 V a (T) = 17 V® (1 - \K^(T)\) 
increases with T for all directions a since the spin correlations tend to zero. For copper 
(Fig. ga)[, 63 V a (T) drops with increasing temperature since the values of the hyperfine 
interactions incidentally are such that the values for fully AFM correlations 63 V^ C (bars on the 
left) are higher than for no correlations 63 V® (bars on the right), and that for both directions 
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a parallel and perpendicular to the planes. The different temperature dependencies of the 
rates measured at the copper and at the oxygen thus find a straightforward explanation. 
While the "basic" relaxation rates are identical for both nuclei, e3 V a (T) has an opposite 
temperature behavior as 17 V a (T) due to the particular values of the hyperfine interaction 
energies, which are atomic properties. 

The temperature dependencies of the rates k U a result from a delicate balance between k V a 
and r e ff(T). From 100 K to 300 K, r e g increases by about 170 %, whereas 63 V C drops by 24 
% but 17 V a increases by 32 %. Some k V a values depend crucially on the precise values of the 
hyperfine energies. Moreover, the rates k T{^ which are directly accessible by experiments 
are linear combinations of the k U a . It is therefore important to study in details the interplay 
between r ej y(T) and k V a (T). 

We illustrate in Appendix [B] the various contributions to 63 V Q (T) and discuss the interplay 
between nearest and further distant AFM spin correlations. 

At this point we just note that in the framework of the MMP model, the Fermi liquid 
contribution which dominates the relaxation rates of the oxygen and yttrium contains no 
correlations, i.e. would correspond to k V® (shown by the high temperature bars in Fig. |HJ). 

It seems now also appropriate to comment on the Korringa relation which has been 
discussed in numerous publications on NMR data analysis in cuprates. The original Korringa 
relation has been derived for an isotropic system. To get a similar relation for layered 
materials would require: (i) temperature independent k Vp and k V-y and a r e ff(T) which is 
proportional to T and to p 2 , the square of the density of states at the Fermi surface, (ii) 
a temperature independent static spin susceptibility in direction a, or (iii) an incidental 
cancellation of the corresponding temperature dependencies in all these quantities which in 
the frame of the present model and in view of Eq. 23 and Fig. 8 is very unlikely. From these 
considerations we conclude that looking for and discussing Korringa relations for NMR data 
in layered cuprates is probably a vain endeavor. 

In Appendix El we show how an apparent linear temperature dependence of the oxygen 
relaxation rate may occur in a limited temperature interval. 
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G. Spin-spin relaxation 



It is worth at this stage to connect the present model for the spin-lattice relaxation in 
cuprates with the theory and experiments of the nuclear spin-spin relaxation rate T^q which 
have been put forward by Slichter and coworker o 32 i 33 i 34 i 35 . T^q measures the strengths of 
the indirect nuclear spin-spin interaction mediated by the non-local static spin susceptibility 
x'if). 

Imai et alM were the first to point out the importance of T^q to obtain information about 
the AFM exchange between the electron spins. They demonstrated that low field NMR 
measurements of 63 T 2 ^ 1 give a strong quantitative constraint on x'{o) i n cuprates. They 
discovered that the staggered susceptibility x'iQ) follows a Curie- Weiss law in YBa 2 Cu 3 7 . 

Since Haase et al^ also presented a derivation of x'(f) m direct space, it is of interest to 
compare Eq. (6) in Ref . S with our expression (j!4b|) for G3 V C (T). The functional form is the 
same but we note that 63 Ti^ requires the knowledge of the correlations at any distance r, 
whereas 63 V C contains only the closest correlations. From the fit of the relaxation rate data 
we got values for K$ x , K X2 and K X3 which, by assumption, could be provided by a static 
spin susceptibility 

X ' a (r) oc cos (Q • f)e- r/Aa . (24) 

We have however no information on x'Jy^) ^ OT 1^1 < 1 and for |r| > 2. Nevertheless, the 
Fourier transform at q = Q gives 

X ' a (Q)=2nXlE X ' a (q = 0) (25) 

where E is an enhancement factor that, in the present approach, cannot be determined 
(expression l25l compares to the isotropic x'{Q) = a £ 2 m the MMP model^). Assuming 
that Eq. (j2IJ) also applies to correlations at arbitrary distances, we evaluate 63 T 2 ^ 1 using 
the values A C (T) obtained from the analysis of the fc T 1 ~ 1 (subsection HVAjl .The resulting 
temperature dependence is shown in Fig. El (solid line) together with the data from Imai 
et al. 3A (circles). F := E\\r = 0)/n% was taken in (|23|) as a parameter to adjust, and 
we calculated it so that our result at T = 300 K corresponds to the data. We found in 
this case F = 8.0 eV -1 . Our prediction of the temperature dependence deviates from the 
measurements, which may indicate that our assumption of an exponential decay of ^'(r) 
over-emphasizes the contributions of high order correlations. If we drop all the correlations 



22 



20 




5- 



TOO 150 200 250 300 

Temperature[K] 

FIG. 9: YBa2Cu307. Comparison with static results (circles) from Imai et alM>. All correlations 
included (full line), only Kq 1} Kqi and K§ x (dashed line) and only NN correlations (dotted line). 

except those coming into the fc T 1 ~ 1 , that is K^, K\ 2 and Kf 3 , we get the dashed line in Fig. 
El with F = 8.6 eV -1 . A better correspondence with the data is obtained if we only keep 
(dotted line), where then F = 11.1 eV _1 . 

On the whole, these results are in agreement with those obtained by Imai et a/ . 33 ' 34 who 
utilized a Gaussian form for the g-dependence of x'io) near q = Q. In this work we have 
used a different parameterization of the AFM correlations (which yields shorter correlation 
lengths), but it seems worthwhile to test further the implications of the present model on 
x'(Q)- An inspection of Eq. 1231 shows that the increase of A^, with decreasing temperature 
determines the increase of staggered susceptibility, provided that Ex'{q = 0) is temperature 
independent. In Fig. EH we plotted the \ 2 a versus 1/T as we obtained them in subsection 
II V Al from the fit on all spin-lattice relaxation data. It is astonishing how well these values 
obey a Curie (or Curie- Weiss) behavior. It seems that upon cooling the spins in YBa2Cus07 
remember their tendency to align antiferromagnetically, but they are prevented from doing 
so as the more energetically-favored superconductivity sets in. This feature has already been 
observed and discussed by Imai et a/ . 33 i 34 . 

An important question at this stage is whether those AFM correlations persist in the su- 
perconducting state. There are strong indications that the correlations develop an extreme 
anisotropy 36 below T c . The question of the persistence of correlations in this state will reveal 
extremely interesting information about the interplay of magnetism and superconductivity 
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FIG. 10: versus 1/T in the direction ab (circles) and c(triangles). Lines are averages, 
which, however, is not the focus of the present work. 

V. UNDERDOPED YBCO COMPOUNDS 

In this section the model is applied to analyze NMR and NQR experiments in 
YBa 2 Cu30 6 6 3 and in YBa 2 Cu 4 8 . We do not have for these compounds the same full 
set of data as for the optimally doped YBa 2 Cu 3 7 , and we therefore restrict our analysis 
to exploring the general trends of the doping and temperature dependencies of the AFM 
correlations and the effective correlation time r e g. We retain the same hyperfme interaction 
energies which have been determined for YBa 2 Cu3C>7 in the analysis of YBa 2 Cus06.63 and 
even for YBa 2 Cu40§. In the latter case, we are able to compare our model predictions with 
high temperature measurements. 

A. YBa 2 Cu 3 6 . 63 

To determine the rates k U a we used the following published spin-lattice relaxation data: 
the copper data from Takigawa et alM (only providing 63 T 1 ~ 1 ), the oxygen data 17 T 1 ~ 1 taken 
from Martindale et alM, and the yttrium data 89 T 1 ~ 1 from Takigawa et al. 27 . They cover 
the temperature range from T c up to room temperature. We are not aware of measurements 
at elevated temperatures. The resulting relaxation rates k U a are shown (symbols) in Fig. 
ITTl The lines in this figure will be explained later. A comparison with those of the opti- 
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FIG. 11: YBa2Cu30g.63- Symbols are experimental data, lines the fit as explained in the text. 
Directions: a (diamonds), b (squares), ab (circles)and c (triangles). 

mally doped material (Fig. ^) reveals that now the temperature dependencies deviate more 
strongly from a linear behavior with a convex (concave) curvature for 63 C/ a (T) ( 17 U a (T)), 
but we do not consider this to be a dramatic contrast. 

The same fitting procedure explained in llV Al was then carried out except that the values 
for the hyperfine interaction energies were kept fixed at the optimized values found for 
YF3a2Cu307 (Table HJ. It is therefore expected that the quality of the fit will be less good 
than for YE^CusCV The result is depicted in Fig. |12(a)[ which shows the temperature 
dependence of the correlation lengths \ a compared with the values for the optimally doped 
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FIG. 12: Correlation lengths and correlations in YBa2Cu30g.63- 



material. As expected, the values for X a are higher (by about a factor 2) in the underdoped 
compound but exhibit a peculiar crossover when the temperature drops below 100 K. In 
this region we also get A c < A a &. Whether this behavior is physically really significant or 
just shows the inadequacy of the postulated model is at the moment open. As was done for 
YBa2Cu307, the correlation lengths can be fitted with a sum of exponential functions. The 
correlations built using the exponential fit are plotted in Fig. |12(b)[ For completeness we 
have plotted in Fig. H3]the calculated k V a . 

The consistency check in analogy with IIV Bl gives 6 different values for k r a (T) which are 
gathered in Fig. Although this time the spread among the different values is considerably 
larger than that for YE^CuaOy (see Fig. |5J) the agreement is still surprisingly good. The 
temperature dependence of the calculated mean values (circles in Fig. [TJJ could not be 
fitted with the function (|21|h However, the same ansatz 



provides a very good fit (as is shown in Fig. IT5|l with the values a = 10 xlO 18 s/K, g= 97 
K, and r 2 = 2.9 x 10~ 15 s. 

The lines in Fig. show the model-calculated U s compared with the experiments, 




(26a) 



but with T\ = aTe 9 ' and r 2 = constant 



(26b) 



26 



1000 = 



-e- 



-e- 



100 



200 

Temperature [K] 



300 



100 



200 

Temperature [K] 



300 



(a) Va 



0.03 



(b) 17 v t 



^ 0.02 



0.01 




100 



200 

Temperature [K] 



300 



(c) 89 V a 



FIG. 13: YBa2Cu306.63- Calculated V a , in the directions a (diamonds), b (squares), ab (circles) 
and c (triangles). 



whereas Fig. ITol depicts the data and the fit in the usual ( k Ti a T) versus T representation. 
The agreement is now only approximate but could be improved by further adjustments 
(in particular the hyperfine field constants) and weighting of the data according to their 
precision. A further confirmation that the fit is however convincing is given in Fig. El where 
we have plotted the ratio 63 T 1 ~ 1 / 17 T 1 ~ 1 . On top of our original data (stars and plusses) and 
our model predictions (full and dotted line), we also show the data (squares) obtained by 
combining the high temperature YE^CusOy measurements of Barrett's et a/.— and Nandor's 
et alM already discussed in II V El At high temperature, we expect the AFM correlations to 
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FIG. 14: r e ff in YBa2Cu30g.63 obtained from a A minimization at the optimized parameters 
(Table IJ). The circles are the calculated average. 17 r a (thin full line), 17 (thin dotted), 17 r c (thin 
dashed), 63 r a b (thick dash-dotted), 89 T a b (thick dash-dot-dot-dashed), 89 r c (thick dot-dash-dash- 
dotted). 
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FIG. 15: T e jj in YBa2Cu3 06.63 obtained from a A minimization at the optimized parameters. The 
circles are the calculated average, and the dashed lines the upper and lower bound. The line is the 
fit as explained in the text. 

become negligible. In this case, the hyperfine fields must be fully uncorrelated and the ratio 
m T^ c l / l7 T{ c l tends to (A 2 ab + AB 2 )/(C 2 a + C 2 ), which is about 8.2. This value is marked by 
a dotted line in Fig. El As can be seen in this figure, the high temperature data confirms 
this prediction, which is a known result for the relaxation by randomly fluctuating hyperfine 
fields. 
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FIG. 16: YBa2Cu306.63- Symbols are experimental data, lines the fit as explained in the text. 
Directions: a (diamonds), b (squares), ab (circles) and c (triangles). 



B. YBa 2 Cu 4 O s 



In stoichiometric YB^CuziOg, NMR and NQR lines are much narrower than in other 
cuprates and allow precise measurements. Raffa et al£L reported high accuracy 63 Cu NQR 
spin-lattice relaxation measurements on 16 and 18 exchanged samples of YBa2Cu408- 
They analyzed their data with the help of the phenomenological relation 



( 63 T lc T) 1 = CT~ a 



1 - tanh 2 (—] 
\2T) 



(27) 
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FIG. 17: ^T^ 1 / 17 T{ c l in YB a2 Cu 3 07 (stars *) and YBa 2 Cu 3 6 . 6 3 (plusses +). Symbols are 
experimental data as cited in the text. The full and dotted lines are the predictions derived from 
the model, and the constant dashed line is the expected value in the absence of AFM correlations. 

which worked reasonably well. To investigate a possible isotope effect of the spin gap pa- 
rameter A they could considerably improve the agreement between data and fit function by 
slightly adjusting the temperature dependence of the function (J2ZJ). 

In terms of our model the temperature dependence of the relaxation rate is given by 

63 T -i = 2 ™ Uah (T) = 2 63 V ab (T) r eff (T). (28) 

To test the quality of the ansatz ()26a|) - (j26b|) for the combined effective correlation time 
we simply take the values for the 63 V a b(T) which we have extracted for the YBa2Cu307 
system. Eq. f!28|) is then used to determine the r e ff(T) from the data of Raffa et a/.—. 
The data (for the 16 sample) and the fit are shown in Fig. |18(a)[ Fig. |18(b)| shows the 
same data with the calculated values extrapolated outside the range (100 K < T < 310 K), 
with the addition of high-temperature data from Curro et alM and Tomeno et al. 39 . The 
temperature dependence of r e g is represented in Fig. [THl The fit gives a = 44 x 10~ 18 s/K, 
g = 195 K and r 2 = 1.9 x 10~ 15 s. Note that these values will change if the 63 V a b(T) values 
appropriate to YBa 2 Cu40g (which are expected to be somewhat larger between 100 K and 
300 K due to enhanced correlations) can be extracted. We guess that our ansatz ()2(iaj) and 
(J26bj) might find a better physical motivation than ()27)1 . Moreover the fit is excellent and 
the extrapolation to higher temperatures is in addition in good agreement with the data 
available. 
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FIG. 18: (a): data from Raffa et al— (crosses) and fit as in the text (solid line), (b): Extrapolation 
of fit to higher temperatures and data from Curro et alm& (triangles right) and Tomeno et al£^ 
(triangles down). 
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FIG. 19: T e ff (circles) is obtained from 63 V^,(T) and the data from Raffa et al. 37 . Dotted line: fit 
as in the text. 

VI. LASCO COMPOUNDS 

A large amount of NMR and NQR data exist also for L^-zSr^CuC^ for various doping 
levels x. Mostly, relaxation rates have been reported for 63 T-^ 1 with only few measurements 
on the oxygens. Haase et alM have been able to determine Kqi(T) (poi(T) in their work) 
from the linewidth of apical, planar O and Cu in La2_ x Sr x Cu04. For the optimally doped 



31 



150 



T, 



'© 



X. 100 





X 







0.1 

Sr content x 



0.2 



50 



-%0 -200 -100 



) 1 
Temperature[K] 



100 



200 



300 



FIG. 20 : 63 Ti c T vs T plot with the data from Ohsugi et a/. 41 for La 2 _ x Sr :r Cu04 . x=0.24 (full 
squares), x=0.2 (empty squares), x=0.15 (full circles), x=0.13 (crosses), x=0.1 (empty circles), 
x=0.075 (full triangles). In the inset, To is plotted versus the strontium concentration. 

compound they find that the correlation decreases with the temperature and is about —0.4 
at room temperature. The correlations are of course expected to be large for low doping. 
In the framework of our model we are not in a position to determine the correlation lengths 
X a b and A c since we do not have enough data. We therefore restrict ourselves in this section 
to a more qualitative discussion of the temperature dependence of 63 T 1 ~ 1 . In particular, we 
can compare the high temperature limit of the effective correlation time r 2 to the relaxation 
mechanism of local magnetic moments in a paramagnet. 

We reproduce in Fig. EH the measurements of Ohsugi et a/.— who present their data in 
a 63 T lc T versus T plot. 

The straight lines that fit the data well at temperatures above 100 K correspond to 



and led, in analogy to the Curie- Weiss law for the susceptibility in insulating antiferromag- 
nets, to the notion of "Curie- Weiss behavior of ( 63 T lc T) " . 

In our model we also expect to find for La 2 _ x Sr :r Cu04 a r e jj as given by (|26aj) and (|26b|) . 
that is 



T lc T = a + aT/T c 



(29) 
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T 3> g. We find indeed in this linear T dependence and can identify 

a ^^VTa ^ T ^~a (31) 

It is astonishing, however, that this linear temperature dependence dominates already at 
T > 100 K. It is probable that like in YB^CusOy, the temperature dependence of m V ab 
is weak. We have shown in Appendix El that due to the incidental values of the hyperfine 
energies for copper, a cancellation of contributions from NN and further distant correlations 
occurs. 

Imai et al.— have studied L^-xSr^CuC^ by NQR up to high temperatures and found 
two striking results. Firstly, at high temperatures (~ 600 K) the 63 T^ rates of all samples 
(with doping levels x= 0, 0.02, 0.04, 0.045, 0.075, and 0.15) were nearly identical although 
the system with x = is an insulator and its magnetic behavior can well be described by a 
two-dimensional Heisenberg antiferromagnet, whereas the samples with x = 0.075 and 0.15 
are conducting and, at low T even superconducting. Secondly, they observed that at high 
temperatures 63 T lc becomes independent of temperature. 

The first striking result of Imai et al—, that is the nearly identical values of m T^~ c for 
x=0 up to £=0.15, is a strong argument for identifying r 2 as a time scale dictated by the 
dynamics of the antiferromagnetism which persists well into the doped and overdoped region. 
This can be illustrated as follows. Moriyap 3 - has calculated T-f 1 for a nucleus of a magnetic 
ion in an insulating antiferromagnet. In the paramagnetic regime he treated the exchange 
interaction within the model of a Gaussian random process with correlation frequency uo e . 
This was adapted by Imai et alr& to a square lattice and appropriate hyperfine interaction 
energies for pure La2Cu04 with the result (equation (4) in Ref. 
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1 



63 T^ = (Al + 4B 2 ) _ (32) 



where cu' e = to e [2A ab B/ (A 2 ab + AB 2 )} 1 ' 2 slightly modifies to e (u' e » 0.9 cu e ). 
In the limit of high temperature, our model simplifies to 

Urn ^ = 2 63 C r 2 = ^ (A 2 ab + AB 2 ) r 2 (33) 



and therefore the relation 
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FIG. 21: La2~ z Sr z Cu04 data from Imai et al. 42 (full squares) and Ohsugi et al4i (for x = 0.15, 
empty squares). 

connects our correlation time r 2 to the correlation frequency of a magnetic moment in the 
paramagnetic regime of an insulating antiferromagnet. It must be stressed, however, that in 
the present work the time r 2 had been introduced as an independent contribution to r e ff(T), 
i.e. r~ff = r^ 1 + t% , which gave an excellent fit to the relaxation rate data in metallic 
YBa 2 Cu307 for 100 K < T < 300 K, with T\ linear in T and r 2 constant. 

To complement these qualitative findings with a quantitative result we present in Fig. |^ 
the data reported by Imai et alM (full squares) and Ohsugi et a/.— (empty squares) for the 
£=0.15 samples together with a fit of the data from the former source, according to (J26bjl - 
()26aj) . The fit, for which we have assumed a constant value for 63 V® b , gives 63 V® b a = 57.1 
(sK) _1 , 63 V® b T 2 = 1574 s _1 and g = 64.1 K. Adopting the hyperfine interaction energies for 
La 2 Cu0 4 reported by Haase et alM, we find a = 106 x 10 -18 s/K and r 2 = 2.92 x 10~ 15 s. 
It is interesting to note that this value of r 2 is very close to those obtained for YBa 2 Cu307 
and YBa 2 Cu30 6 . 63 . 

VII. SUMMARY AND CONCLUSIONS 

We have compiled a complete set of relaxation rates data for optimally doped YBa 2 Cu307 
from T c up to room temperature and extracted k U a , the contributions to the rates from 
the fluctuations along the individual crystallographic axes. The result of this simple data 
transformation is that from this point of view there is no striking difference between the 
copper and the oxygen relaxations. This suggested the decoupling k U a (T) = k V a {T)r e ff(T), 
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where the temperature dependence of k V a is given by the static AFM correlations (T) 
for the oxygen and Kq X {T), K° 2 (T) and K° 3 (T) for the Cu. A numerical fit to the data on 
•^-independent ratios of relaxation rates yielded the correlation lengths in function of the 
temperature as well as refined hyperfine field constants. The validity of our approach is 
confirmed by the fact that the seven h r a (T) are well grouped. The average, assimilated to 
the effective correlation time r e ff, is extremely well modeled by r~^(T) = r 1 _1 (T) + r 2 _1 with 
r 2 constant and t±(T) linear in T. This led to the surprising result that in YBa 2 Cu 3 07 the 
"basic" relaxation of all three nuclei under consideration is dominated at low temperature 
by scattering processes of fermionic excitations. The extrapolation of the model predictions 
to higher temperatures is in very good agreement with the measurements. A similar but 
slightly reduced analysis was conducted on YBa 2 Cu 3 6 . 6 3, leading as expected to higher 
values for the correlations. In the case of YE^CuztOs we had an even more reduced set 
of data, however some very precise data for Cu are available at high temperature. The 
analysis of the effective correlation time T e g{T) in both underdoped compounds revealed 
that T~ff(T) = rf^T) + t 2 1 fits the result again very well, provided that t\(T) is modified 
by a gap function at lower temperature. From the analysis of the La^zSr^CuC^ series we 
could connect r 2 with relaxation due to AFM spin fluctuations in a paramagnetic state. 

In conclusion, we have shown that in the model of fluctuating fields the AFM spin-spin 
correlations Kq^T), K^ 2 (T) and K^(T) determine the degree of coherency. In our fit 
to the data we found that the in-plane correlations are about —0.4 for YBa 2 Cu 3 7 and 
—0.65 for YBa 2 Cu408 at 100 K. The copper and yttrium relaxation rates depend on all 
three correlations Kq^T), K^ 2 {T) and K" 3 (T), whereas 17 T 1 ~ 1 involves only Kq ± . Moreover, 
the contributions of the three correlations at the copper partially cancel out due to the 
particular values of the hyperfine field constants. It would be therefore extremely useful 
to collect more measurements on oxygen, in all crystallographic directions and in all 
substances. 

We also note that whereas the particular choice of spatial exponential decay of the 
correlations fixes the details of the correlation lengths and hyperfine field constants, the 
same general conclusions of our model would be drawn if another form of spatial decay 
would be applied. Another note is that the question of the degree of coherency should be 
addressed at an even lower level than we do, since one should also consider how the various 
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contribution to the on-site hyperfine fields are added together. For simplicity however we 
postpone this discussion to a future publication. 

The essential asset of the model presented in this paper is that a same (fast) fluctuation 
mechanism of the spin liquid can be identified in the relaxation of all the nuclei (copper, 
oxygen and yttrium), and this mechanism is characterized by a temperature-dependent 
effective correlation time r e ff(T). Moreover, the parameterization T e g~ x = r i _1 + r 2^ 1 has 
a wide range of validity. At low temperature, the effective correlation time r e g is just 
Ti, which is linear in YBa2Cu3C>7 but is modified by a gap function in the underdoped 
compounds. In YE^CuaOy the term t% could be linked to the nuclear relaxation by charge 
carriers in metals. At high temperature, the effective correlation time r e ff goes over to 
the constant T2, which could be connected to the correlation time in antiferromagnets. 
This seems to indicate that at high temperature, the nuclei are probably influenced by a 
system of local moments not very different from the paramagnetic phase of the underdoped 
parent compounds. Lowering the temperature, a smooth crossover to an itinerant system 
is observed. It looks like these features which have been observed by Imai et alM in 
La2Cu04 are also present in the YBaCuO system although the crossover happens at higher 
temperatures. At low temperature, even in the optimally doped YBa2Cus07, a seemingly 
local pairing of moments occurs that is reflected in the Curie- Weiss behavior of x'(Q)- 
Further analysis, however, is required to get more quantitative results about the doping 
dependence of all these phenomena. It should be mentioned that also heavy fermions and 
mixed valent systems exhibit a crossover from localized moments to coherent behavior with 
decreasing temperature, as has been pointed out recently 45 . 

We also note that the incommensuration or discommensuration of the AFM ordering 
observed by neutron scattering measurements does not invalidate the model presented 
here 5 ^. Moreover, the model ()26a)l and (J26b)l has a wide range of application, extending 
also to electron-doped superconductors. In particular, the NMR measurements taken by 
Imai et alM on Sro.gLao.iCuC^ can also be fitted within our model. 

Another important conclusion is that having identified a similar "basic" nuclear 
relaxation rate for all nuclei, the different temperature behavior of the observed rate of 
the copper and oxygen can be put down to the particular values of the hyperfine field 
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constants. A detailed discussion of the functions k V a (T) showed the crucial interplay of 
nearest neighbors and higher AFM correlations and revealed the complications one meets in 
trying to find Korringa relations in layered cuprates. A case of practical interest is that the 
on-site and transferred hyperfine energies for Cu are such that 63 V a b(T) (which determines 
63 T 1 ~ 1 ) changes only slightly since the contributions from the NN and the third neighbor 
correlations nearly cancel each other out. 

We proposed in this paper an alternative way of looking at the spin-relaxation data 
which we believe could help identifying the underlying physical processes in the cuprates. 
The analysis of the selected set of data in cuprates shows that the present model reveals 
the wealth of information that can be obtained from NMR and NQR measurements. While 
the success of the parameterization t~^(T) = rf *(T) + r^ 1 is astonishingly good and 
universal, we could only suggest likely explanations about its significance, and a conclusive 
interpretation is still open to speculations. 
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APPENDIX A: SCALED CONTRIBUTIONS TO THE RELAXATION RATES 

The quantities k U a extracted from the experiments can be gathered into one plot, each 
on a different scale. This is equivalent to applying the affine transformation 

k U a = k f a k U a + k g a (Al) 
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(a) YBa 2 Cu 3 7 (b) YBa 2 Cu 3 6 . 63 

FIG. 22: The k U a after afflne transformations of the data. 17 U a (thin full line), 17 Ut> (thin dotted), 
17 U c (thin dashed), 63 U ab (thick dash-dotted), 63 I7 C (thick dashed), m U ab (thick dash-dot-dot- 
dashed), 89 U c (thick dot-dash-dash-dotted). 



with constants k f a and k g a . The result for YBa2Cu3C>7 is shown in Fig. |22(a)[ where the 
vertical axis is arbitrary. We note that the k U a in YBa2Cu3C>7 have the same tempera- 
ture dependence for most of the temperature range. The larger deviations occur at higher 
temperature for 63 U a b (thick dash-dots) and 17 Ub (thin dots). The same operations can be 
done for YBa2Cu306.63, Fig. |22(b)| Again all k U a fall on the same curve, apart 63 U a b which 
deviates significantly. 

A similar picture is drawn upon dividing the rates k U a (T) by the constants k V® that are 
calculated from fl"5|) and whose values are gathered in Table HU k U a (T) := k U a (T)/ k V® ■ 
The scaled relaxation rates k U a (T) are shown in Fig. |2Hlfor YBa2Cu3C>7 and YBa2Cu30g.63- 
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(a) YBa 2 Cu 3 07 



(b) YBa 2 Cu 3 6 .63 



FIG. 23: k U a (T) = k U a (T)/ k V°. 17 U a (thin full line), 17 U b (thin dotted), 17 U C (thin dashed), 
63 U ab (thick dash-dotted), e3 U c (thick dashed), m U ab (thick dash-dot-dot-dashed), 89 U c (thick dot- 
dash-dash-dotted) . 



APPENDIX B: DETAILED INVESTIGATION OF k V a (T) AND KORRINGA- 
LIKE BEHAVIOR 

We plotted in Fig. |2Uthe individual contributions to 63 V a (T), which are defined as 

63y a(T) = ^A a BK^{T)\ 
63 V*(T) = i[8B 2 ^(T)] 

m V 3 {T) = ^- 2 [4B 2 K? 3 (T)} (Bl) 

63 V^, 63 and 63 V 3 are the terms proportional respectively to the AFM spin correlation 
functions Kq X , K" 2 and K± 3 . We also have the constant term 63 V®, the temperature 
independent contribution given in Eq. (|15a|) and Table HU Since A ab > 0, A c < 0, 
< and K" 2 ,K^ 3 > 0, all individual contributions 63 are positive except 63 V^ b . 
We see in Fig. |24(a)| that 63 V ab (large circles) is made up almost entirely of 63 V^ + 63 V a 2 b 
(small circles), whereas 63 V c in Fig. |24(b)| (large triangles) is dominated by 63 V C ° + 63 V^ 
(small triangles). In other words, in the ab-direction 63 V* b and 63 V 3 b , having opposite 
signs, almost cancel each other out. For example near T c where the correlations are high, 
I 63 Kb/ 63 Kb I = 2 ^l-^oil/l-^oil 2 ~ 1- m tne c-direction all contributions have the same 
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FIG. 24: YB a2 Cu 3 07. Full m V a (large symbols), m V® (dotted line), 63 V£ (full line), m y£ (dashed 
line), 63 V* (dash-dotted line). 63 V$ + 63 V* b (small circles), 63 V C ° + (small triangles). 
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FIG. 25: YBa2Cu307. Decomposition of (1 + K$)T e jj (full line): t\ (dashed line), r e jj — t\ 
(dash-dotted line) and K^T e g (dotted line). 



sign but due to the high value of |A C | the lower order correlation plays the major role. 
Note that the same qualitative behavior is expected irrespective of the choice of correlation 
dependence (Eq. H7)l . 

The seemingly Korringa-like behavior of ( 17 Ti a T) 1 is investigated in Fig. 123 There we 
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decomposed 17 T lc 1 oc (1 — |-K"oil) r ej8 F (full li ne ) i n t° 

(1 - |<|)r eJf = n + (r eff - n) - r ef |<| (B2) 

From the figure 051 we see that the contribution T e g\KQ\\ (dotted line) varies little over a 
large range of temperature. The dashed line is the truly linear contribution T\. The actual 
temperature dependence of ( 17 Ti a T) 1 is thus given by this linear contribution, minus that 
of (r e ff — Ti) (dash-dotted line). The total result has however the appearance of linearity. 
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